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ABSTRACT 



(N 
> 

On 

f^* ) Context. Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal 
structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material 
\Q in the immediate vicinity of the central black holes in AGN. 

, Aims. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well-known 
*0 quasar 3C 273. 

' Methods. We present here the results from a 3D relativistic hydrodynamics numerical simulation based on the parameters 
C"| given for the jet by Lobanov & Zensus (2001), and one in which the effects of jet precession and the injection of discrete 
O |' components have been taken into account. We compare the model with the structures observed in 3C 273 using very long 
baseline interferometry and constrain the basic properties of the flow. 

Results. We find growing perturbation modes in the simulation with similar wavelengths to those observed, but with a different 
set of wave speeds and mode identification. If the observed longest helical structure is produced by the precession of the flow, 
, longer precession periods should be expected. 
• • i Conclusions. Our results show that some of the observed structures could be explained by growing Kelvin-Helmholtz instabilities 
^ ' in a slow moving region of the jet. However, we point towards possible errors in the mode identification that show the need of more 
complete linear analysis in order to interpret the observations. We conclude that, with the given viewing angle, superluminal 
» I ■ components and jet precession cannot explain the observed structures. 
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1. Introduction 

The structure and kinematics of parsec-scale outflows is 
typically explained in terms of shocks (Marscher [1980 
Marscher & Gear IT5g51 Gomez et al. IT5531 IT551 Gomez 
et al. I1994[) and Kelvin-Helmholtz (K-H) instabilities 
(HardeeUnmUnHlinHa Hardee et al. 115571 HardeeEOOO, 
2003, Hardee et al. l2005fl developing in a relativistic fluid. 
Relativistic shocks may dominate the jet dynamics and 
emission at small scales, but are likely to dissipate at 
distances larger than ~ 10 pc (Lobanov & Zensus H999[) 
due to the interaction with the slower flow. On interme- 
diate scales (~ 10-100 pc) shocks and plasma instabili- 
ties may play equally important roles in jets (Lobanov & 
Roland 2001). Distributions of the synchrotron turnover 
frequency obtained for 3C 273 (Lobanov et al. 1997) and 
3C345 (Lobanov 1998) indicate that both shocks and in- 
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stabilities are present on these scales, while larger scales 
are most likely dominated by plasma instabilities alone. 

Recent studies by Hardee (J2000) have shown that 
K-H instability may produce complex, three-dimensional 
ribbon-like and thread-like patterns inside a relativistic 
jet. In these ribbons and threads, a substantial increase 
of particle pressure and radio emissivity can be expected. 
This model has been successfully applied to the jet in 3C 
120 (Hardee l250H Hardee et al. 1255511 . The threaded struc- 
ture forming a double helix has been detected in a space 
VLBI image of 3C 273 made at 5 GHz (Lobanov et al. 
200(| . It was explained in terms of K-H instability devel- 
oping in a relativistic flow with a modest Lorentz factor 
7 = 2.1 and a relativistic Mach number M = 3.5 (Lobanov 
& Zensus 2001, hereafter LZ01). The analytical approach 
used in LZ01 is based on linear perturbation analysis of 
a K-H instability developed by Hardee l(T557l I2555|) . A 
similar approach applied to kiloparsec-scale jet in M87 al- 
lowed for accurate determination of physical parameters 
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and modelling of radio emission to be made (Lobanov et 
al. 1200311 . 

However, results from numerical simulations of rela- 
tivistic flows indicate that, after the linear regime of in- 
stability growth, the jets can be easily disrupted (Perucho 
et al. 1200 4b). In addition to this, the bulk Lorentz factor 
7 = 2.1 derived in LZ01 is below the values required to 
explain the apparent speeds of /3 apP ~ 5 — 8 c of enhanced 
emission features observed in 3C273. LZ01 suggest that 
the K-H instability is developing in a slower, underlying 
flow, and the fast components are most likely faster shock 
waves produced in the jet by the periodic ejections associ- 
ated with the nuclear flares. The presence of such shocks 
may disrupt the linear growth of the K-H instability as 
well, and it is not clear whether the linear stability anal- 
ysis can be still applied in the presence of these kind of 
non-linear effects in the flow. 

In view of these concerns, it is important to confront 
the results of LZ01 with numerical simulations, and at- 
tempt to address several fundamental issues about the 
stability and propagation of relativistic flows similar to 
the one observed in 3C273. Numerical simulations can 
be used to verify whether the linear theory can be ap- 
plied for explaining self-consistently the morphology and 
kinematics of parsec-scale flows, and whether these flows 
preserve fingerprints of linear modes even when the non- 
linear regime has developed or non-linear features, such 
as fast components, appear. Numerical simulations pro- 
vide a means to address these problems by following the 
transition from linear to non-linear regimes of instability 
development (Perucho et al. 12004a, 2004b|. The ultimate 
goal of this work is to probe the advantages and limi- 
tations of the combination of different approaches (linear 
theory, numerical simulations and observations) to studies 
of parsec-scale jets. 

3C 273 is the second quasar discovered (Hazard et al. 
1963), and the first one for which the emission lines were 
identified with red-shifted hydrogen lines (Schmidt |T963). 
In the same work, Schmidt (1963) also pointed out the 
presence of a jet-like structure in this object. During the 
last four decades, the active nucleus and the relativis- 
tic outflow in 3C273 have been studied in great detail 
(Courvoisier 1998). The parsec-scale radio jet in 3C273 
has been monitored for almost three decades (Pearson et 
al. rrWn Unwin et al . [TTiSSI [TTiSOl Zensus et al . [T3EBI ITOOOI 
Davis et al. 119911 Abraham et al. 119961 Krichbaum et al. 
120001 Lobanov et al. 120001 Asada et al. l2002|l . The emis- 
sion associated with the relativistic outflow on kiloparsec 
scales has been probed extensively in the radio (Conway 
et al. [TOSH H333)l . near infrared (Neumann et al. HMTl 
Hutchings et al. I2001J) . optical (Thompson et al. H3331 
Jester 120011 Jester et al. I2001|) and X-ray (Roser et al. 
120001 Marshall et al. I200T1 Sambruna et al. I200T) wave- 
bands. 

The relativistic jet observed in the quasar 3C273 is 
one-sided, with no signs of emission on the counter-jet 
side at dynamic ranges of up to 16,000:1 (Unwin et al. 
I1985|) . This is evidence for strong relativistic boosting in 



an intrinsically double-sided outflow powered by an ac- 
cretion disk around the central black hole (Begelman et 
al. HSH1) . The mass of the central black hole in 3C273 
is estimated to be M b h = 5.5to;§ x 10 8 Mq (Kaspi et al. 
2000). The enhanced emission features (jet components) 
identified in the jet on scales of up to ^20 milliarcseconds 
(mas) are moving at apparent speeds exceeding the speed 
of light by factors of 5-8 (Abraham et al. I1996[) . Plausible 
ranges of the Lorentz factor 7 w 5-10 and viewing angles 
9j e t « 10°— 15° have been inferred from these measure- 
ments. 

Ejections of new components into the jet occur roughly 
once every year (Krichbaum et al. I2000J) . and they are 
likely to be related to weak optical flares observed with 
a similar periodicity (Belokon 1981). The position angle 
at which the components are ejected shows regular varia- 
tions with a likely period of about 13-15 years (Abraham 
et al. 119961 Abraham & Romero 1999), correlated with 
the long-term variability observed in 3C 273 in the optical 
(Babadzhanyants & Belokon 1993) and radio (Turler et al. 
1999| bands. Abraham & Romero (1993 have suggested 
that this periodicity may reflect changes of the jet axis 
induced by the relativistic precession of the inner part of 
the accretion disk. 

Results from the linear analysis and numerical mod- 
elling are presented, compared and discussed in Sect.[2Jin 
connection to explaining the observed properties of parsec- 
scale outflow in 3C273. Main results of the work are dis- 
cussed in Sect. 01 

Throughout the paper, we adopt the flat ACDM 
Cosmology with the Hubble constant Hq = 71 h km 
S _1 Mpc _1 , where h is a constant with a likely value of 

1, and matter density £Im — 0.27. The positive definition 
of spectral index, S cx v a is used. For 3C273 (z — 0.157, 
Strauss et al. 1992}, the adopted cosmological parameters 
correspond to the luminosity distance Z?l = 0.7ft. -1 Gpc. 
The respective linear scale is 2.69/i _1 pemas -1 , and a 
proper motion of lmas/yr corresponds to an apparent 
speed of lO.l/i -1 c. 

2. Numerical simulations 

In this section, we present numerical simulations per- 
formed with the aim to provide a counterpart to the an- 
alytical modelling of the jet structure in 3C 273 made 
in LZ01. To give a better account of the connection be- 
tween the analytical and numerical approaches, basic re- 
sults from the linear model of LZ01 are briefly summarized 
below. 

2.1. Results from the linear analysis 

Linear perturbation analysis of Kelvin-Hclmholtz insta- 
bility (cf. Hardee 1T0871 120001 and Hardee et al. lT097| was 
applied in LZ01 to explaining the internal structure of the 
jet in 3C273. The locations of two thread-like features 
identified inside the jet (the features PI and P2, in the 
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nomenclature of LZ01) were approximated by combina- 
tions of several oscillatory modes. These modes were iden- 
tified with different modes of Kelvin- Helmholtz instability. 
The parameters of these modes are given in Table Q] The 
characteristic wavelengths of different instability modes 
can be related to the jet speed ((3j), Mach number (Mj) 
and the density ratio (77 = Pj/p a ) between the jet and the 
ambient medium. These wavelengths are obtained by ap- 
proximating the dispersion relation in the small frequency 
limit (to — > 0) and large frequency limit (u> 3>). Both are 
obtained for highly supersonic jets (Mj 3> 1). The longest 
unstable wavelength, X l nm (obtained in the low-frequency 
limit), and the resonant wavelength, A* m (obtained in the 
high frequency limit), are given by fHardee ll987[l : 

4 7j ^(AfJ- 1)1/2 

n + 2m-l/2 ' { ' 

a* 2lT 

Anm /3., a //3j (n + m+ 1/2) 

(M2-/3])V2 +7i (M|-^)V2' [> 

where f3 s , a is the sound speed of the external medium in 
units of the speed of light, jj is the Lorentz factor of the 
jet, M a = Pj/Ps,a, Rj is the jet radius, and n (the az- 
imuthal number) and m give the number of nodes around 
the jet surface and the number of nodes between the axis 
and the surface, respectively. The first equation gives the 
longest unstable wavelength for a body mode (m > 0) 
and the zero frequency limit for a surface mode (m = 0), 
as the latter do not show the long wavelength cut, and 
the second stands for the most unstable wavelength of a 
given mode. A total of five wavelengths were identified 
from fits to the double ridge line found in the observa- 
tions presented in LZ01. Physical parameters of the jet 
were obtained in LZ01 by relating these wavelengths to 
the wavelengths of the oscillatory modes from the fit to the 
internal structure of the jet. The jet parameters obtained 
are given in Table |21 Four of the observed wavelengths 
were associated with resonant wavelengths of helical and 
elliptical modes and the longest wavelength (18 mas) was 
associated with the helical surface mode driven externally 
at approximately twice the resonant wavelength. 

It should be noted that the Lorentz factor 7 = 2.1 
derived for the jet is below the values given by other au- 
thors in order to explain superluminal motions observed 
in 3C 273. This may result from Kelvin-Helmholtz insta- 
bility developing in an underlying, slower flow, and not 
in the flow that contains ballistic, superluminal features 
(LZ01). 

The results from the linear stability analysis are com- 
pared to the numerical solutions of the stability problem 
obtained from the individual simulation runs described 
below. 



Table 1. Identified wavelengths, modes and their ampli- 
tudes from observations (LZ01). H stands for helical, E 
for elliptical modes, and subscripts refer to surface (s, fun- 
damental) or body modes (b, reflection). The latter are 
followed by the index identifying the exact body mode.* 
stands for identified resonant modes. 



A [mas] 


Amplitude [mas] 


Mode 


PI P2 


PI P2 




18 


1.5 


H e 


12 


1.4 


El 


3.9 4.1 


2.2 1.5 




3.8 


1.2 


E'ti 


1.9 


0.25 


tlbi 



2.2. General properties of the numerical models 

Numerical simulations were performed using a three- 
dimensional finite-difference code based on a high- 
resolution shock-capturing scheme which solves the equa- 
tions of relativistic hydrodynamics written in conservation 
form. This code is an upgrade to 3D of the code described 
in Marti et al. 1)1997(1 and shares many features with the 
3D code GENESIS (Aloy et al. H333} . It was parallelised 
using OMP directives. Simulations were performed in an 
SGI Altix 3000. 

In the numerical models, the initial properties of a 
stationary flow in pressure equilibrium with the external 
medium are set according to the results of the linear mod- 
elling (see Table 0). It has to be noted, however, that no 
opening angle has been taken into account. Perturbations 
are applied at the inlet. Boundary conditions are: 1) in- 
jection at the inlet (with the parameters given by LZ01), 
and 2) outflow at the side boundaries and at the axial end 
of the grid. An extended grid with a decreasing resolu- 
tion is added on each side of the main grid and at its axial 
end, in order to avoid spurious numerical reflections of the 
solution at the main grid boundaries. 

In order to achieve a steady initial model, we add a 
smooth transition (i.e., shear layer) between the jet and 
the ambient medium, of the form: 

p(r)=p a -(p a - Pj )/cosh(r) m , (3) 

v z (r)=v z , /cosh(r) m , (4) 

where p stands for rest mass density and v z for axial ve- 
locity (v z ,o is the value at the axis corresponding to the 
Lorentz factor -fj = 2.1), subscripts a and j correspond to 
ambient medium an jet, respectively, and r is the radial 
coordinate. The smaller the resolution, the smaller the ex- 
ponent m has to be in order to reduce the numerical noise 
below the amplitudes of the perturbations. 

We performed two simulations to investigate the gen- 
eral development of a K-H instability in the flow and to 
analyse the effect of fast components and the jet pre- 
cession. In the first simulation (3C273-A) we perturb a 
stationary flow so as to observe which modes and wave- 
lengths dominate the jet structure. In the second Simula- 
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Table 2. Jet parameters from the fit. jj is Lorentz factor, Mj r is the relativistic Mach number, r\ is the jet-to-ambient 
rest mass density ratio, 4>j is jet half opening angle, 6j is jet viewing angle, c s j a are sound speeds, and / is the projected 
linear scale. 



li 


M 3,r 


V 






% [°] 


c 3 ,j [c] c S)a [c] I [pc/mas] 


2.1 


3.5 


0.023 


0.8 


1.5 


15 


0.53 0.08 2.43 



tion (3C273-B) we try to see if similar instability struc- 
tures can be generated by the precession of the jet and 
periodic injections of faster components into the flow. In 
both simulations, the parameters of the steady flow are 
those from LZ01 (see Table |2Jl. We use the perfect gas 
equation of state with the adiabatic exponent Y = 4/3. 

2.3. Simulation 3C273-A 
2.3.1. Initial setup 

In this simulation we introduce perturbations at frequen- 
cies calculated such that they are expected to reproduce 
the observed wavelengths in the jet structure if these are 
propagating at the given speed in LZ01. The grid extends 
over 844 cells in the axial direction and 128 cells in both 
lateral directions (including the extended grid). We use 
a resolution of 16cells/Rj, with Rj the radius of the jet, 
in the transversal direction and 4cells/Rj in the direc- 
tion of the flow. Simulation lasted for a time 1097 Rj/c 
(i.e., about 2852 yrs when scaled to source units; see next 
paragraph), and it used ~ 11 Gb of RAM memory and 
8 processors during around 30 days in a SGI Altix 3000 
computer. 

Assuming an angle to the line of sight of 15° and red- 
shift z = 0.158 (1 mas = 2.43 pc) the observed jet is 169 pc 
long. Considering the jet radius given in LZ01 (0.8 pc), 
the numerical grid extends over 211 7?-, (axial) times 8 Rj 
times 8Rj (transversal), i.e., 169 pc x 6.4pc x 6.4 pc. This 
allows us to accommodate all relevant relativistic and sub- 
relativistic structures which could give rise to the wave- 
lengths observed in the patterns PI and P2 from LZ01. 
A shear layer of 2 Rj width (to = 2) in Eqs. is 
included in the initial rest mass density and axial velocity 
profiles to keep numerical stability of the initial jet. To 
avoid reflection of the numerical noise from the main grid 
boundaries, an extended grid is introduced, with a cell 
size increasing progressively by 20%, as we do not need 
fine cells in this extended region because we are mainly in- 
terested in the linear regime, which does not involve large 
radial distortions of the jet. The extended grid has 24 cells 
in the radial directions, reaching up to 36 Rj on each side 
of the jet, and 168 cells axially, reaching up to 316.5 Rj. 
Elliptical and helical modes are induced at the inlet using 
the following expression: 

P' = A ° 2 cos(ujt + n6) sin 2 (tit) , (5) 
cosh r 

where A (= 10 -4 ) is the initial amplitude, r is the radial 
coordinate, uj is the frequency of the mode, n — 1 for 
helical modes and n = 2 for elliptical ones, 8 is the polar 



angle in cylindrical coordinates, and sin 2 (7rr) is used in 
order to give an initial transversal structure to the modes. 
The evolution of perturbations and their coupling to K- 
H modes have been shown to be independent from this 
transversal structure (Perucho et al. !2005f) . In Fig. 2] we 
show axial and transversal cuts for a three-dimensional jet 
with periodic boundary conditions in the axial directions, 
like those used in Perucho et al. I|2005|l . The axial structure 
shown in the axial cuts is added to Eq. JSJ as a term 
k z z in the cosinus function. Typical structures induced 
by Eq. JSJ) in such a generic jet are those shown in the 
transversal cuts. The sum of all the input modes gives 
the total perturbation. The simulation has to reproduce 
the resonant wavelengths of the basic modes identified in 
Table |2 (2, 4 and 12 mas). It should be noted that the 
helical surface mode at Ah s = 18 mas is driven externally 
and therefore cannot be reproduced in this simulation. 

Frequencies of the excited modes, both helical and el- 
liptical, are introduced in Eq. [3] These frequencies are 
computed from the observed wavelengths in LZ01, A obs , 
corrected for projection effects and relativistic motion of 
the wave (time delay), with velocity v w (— 0.23c). We use 
uj = 2irv w /\ theor (see Table d, where 

x theor = X° bS (l-Vw/c COsOj) 

sin 0j 

where 9j is the angle to the line of sight. We observe that, 
when computing the numerical solutions for the stability 
problem for the proposed jet parameters, the wavelengths 
and speeds obtained for the relevant modes (see Table |2J 
are different to those given in LZ01. This fact may be 
caused by errors in the angle to the line of sight or in the 
derived parameters in LZ01, including the wave speeds 
which are used in the transformation. 

2.3.2. Results 

To compare the results of the simulation with the numer- 
ical solution of the stability problem, we have solved the 
dispersion relation equation (e.g., Hardee 2000) for the pa- 
rameters given in Tabled Fig. [21 shows the solutions for 
the pinching, helical and elliptical surface, first and second 
body modes. Table0]shows the characteristic wavelengths 
of the relevant modes. 

Fig- El presents axial cuts made at two different times 
of the simulation. In Fig. [3] we observe a A ~ 20 — 25Rj 
helical structure in the upper plot that could be associ- 
ated to u) 2 (see Table ISJ) . Fig. 0] shows several transversal 
cuts of the jet illustrating competition between the heli- 
cal and elliptical modes. We can see how excited modes 
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Table 3. Correspondence between the observed wavelength and the perturbations used for the simulation 3C273-A. 
X obs is the observed, projected wavelength, \ theor is the intrinsic wavelength, and ui is the derived frequency for the 
given intrinsic wavelength and wave speed (see text). Fourth to seventh columns give the wavelengths and wave speeds 
derived from the dispersion relation solution, marked in Fig. [3 as w±, u>2 and 0J3. 



\° bs [mas] 




u [c/Rj] 


Pinch bl 


Helical s 


Elliptical s 


Elliptical bl 








A; v w 


A; v u 


A; v u 


A; v w 


12 (E*) 


110.0 


0.013 (wi) 




38.5; 0.081 


18.5; 0.077 




4(^*i,^i) 


37.4 


0.039 (u 2 ) 




16.0; 0.10 


6.1; 0.078 




2 (h; 2 ) 


18.7 


0.078 (u 3 ) 


7.0; 0.088 


12.1; 0.15 


5.4; 0.14 


3.2; 0.08 



Model - BDS.pird- Model = B05,h«Kcal 

Time - 0. R/c Time = 0. R,/c 

Log. Pressure Log. Pressure 




Fig. 1. Structure of perturbations (axial and transversal cuts) in a generic 3D jet, as described by Eq. 10. Top left: 
pinching mode (n = 0). Top right: helical mode (n — 1). Bottom left: elliptical mode (n — 2). Bottom right: triangular 
mode (n = 3). 



Table 4. Solutions of the stability problem in the fastest 
growing or resonant frequencies (indicated with an aster- 
isk, column one lists the frequencies, column two the wave- 
lengths and column three the wave speeds) and longest 
unstable wavelength (fourth column, A;) 



Mode 


u>* (c/Rj) 


A* m 


vZ (c) 


A* (Rj) 


Pinch bl 


0.46 


3.5 


0.26 


7.5 


Pinch b2 


0.86 


1.8 


0.25 


3.1 
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0.24 


7.6 
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Helical b2 
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0.25 
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Elliptic s 


0.19 


5.3 


0.16 




Elliptic bl 


0.85 


1.9 


0.26 


3.2 


Elliptic b2 


1.27 


1.3 


0.25 
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dominate at different positions and times in the jet. It 
is remarkable that elliptical structures show up close to 
the injection point, while helical modes, develop in the 
jet farther downstream. This agrees with the conclusions 
presented in LZ01. Nevertheless, we have not been able 
to clearly identify the elliptical mode in the longitudi- 
nal cuts. Fig. [S] shows longitudinal cuts of pressure per- 
turbation (defined as the difference between the value of 



the pressure in a cell and the initial equilibrium pressure, 
P — Pq, with _P ~ 0.03 p Q c 2 ) at different jet radii, from 
which the dominant wavelengths could be identified in the 
simulated jet. We identify a A ~ 40 — 50 Rj structure at 
z < 60 Rj which we interpret as due to beating between 
two wavelengths of the first body helical mode at wave- 
lengths 4.5 Rj and 4Rj, like that derived from the fits to 
the observations by LZ01, and given in Tabled (3-9 and 
4.1 mas). From plots of pressure perturbation at different 
radii (Fig. |5J) we conclude that the radial structure of this 
beat can only be produced by the helical first body mode, 
as the fluctuations are stronger at Rj/2. The beating could 
also be produced by the elliptical surface mode, but the 
fact that pressure fluctuations are smaller at the jet sur- 
face rules out this possibility (see, e.g., Hardee [2000 ). At 
larger distance (z > 70 Rj) we have identified a large am- 
plitude helical ~ 25 i?, wavelength and a shorter 2.5 Rj 
wavelength superposed on the former one. In an axial cut 
of the pressure perturbation close to the axis we have iden- 
tified an elliptical fluctuation with wavelength ~ 3.5 Rj at 
z ~ 70 — 95 Rj and a helical one with wavelength 2.5 Rj at 
z > 100 Rj. All the modes that are reported in this para- 
graph are pointed in Fig. [5] with arrows indicating their 
wavelengths and wave speeds. We should keep in mind 
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Fig. 2. Solution of the stability problem for the parameters given in TableEJ In the upper plots, real (dashed lines) and 
imaginary (dotted lines) parts of the wavenumber are given as a function of frequency. In the lower plots, corresponding 
wavelengths and wave speeds are shown. Frequencies col, lu2 and w3 correspond to those given in Table [21 for the 
simulation 3C273-A. Note that wavelengths for the elliptical modes in the plot are for a 180° rotation of the ellipse 
and needs to be doubled for a complete rotation. This is because such a 180° rotation generates a repetition of 
the structure in elliptical modes, therefore giving the impression of a complete wavelength. Arrows connecting the 
wavelength and wave speed plots indicate identified modes in the simulation 3C273-A (see text). 



that the stability problem has been solved for the vortex 
sheet case and the jet in the simulations has a thick shear 
layer. This fact can introduce inaccuracies in the detection 
of modes in the simulation. 

We have used all the identified structures in order 
to produce theoretical cuts of the pressure perturbation. 
Results are shown in Fig. H3 Note that here we have plot- 
ted P/Po, whereas in Fig. we plot P — Pq. The similarity 
to that obtained from the simulation is remarkable taking 
into account the fact that the presence of a shear layer, 
which may modify the stability problem solution (Perucho 
et al. l2005[l . has not been considered in the interpretation 
of the results. 

When comparing the results given in the previous 
paragraphs with those from LZ01, we find differences in 
the typical wave speeds for the modes observed in the 
simulation (~ 0.1c), obtained from the solution to the 
stability problem (Fig. |2J, compared to those given in 
LZ01 from the linear approximations (0.23c). This could 
be caused by the uncertainties of the linear approxima- 
tions, Eqs. and @, used in LZ01. In deriving ap- 
proximations, a large classical Mach number is assumed 
(Mj = Vj/c s j 3> 1), but this is not generally the case for 
hot jets, for which c s j ~ 0.57, and thus, Mj < 1.75. To 
investigate the uncertainties introduced by this fact, we 
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Fig. 5. Longitudinal cuts of pressure perturbation at dif- 
ferent radii (Rj/8 top left, 3Rj/8 top right, 5Rj/8 bottom 
left, 7Rj/8 bottom right) and in symmetric positions with 
respect to the jet axis at t = 250 Rj/c. 

have used numerical solutions of the dispersion relation 
for different cases and compared them with the results 
of the approximations. Results show that the errors in 
the determination of characteristic wavelengths with lin- 
ear approximations can reach a factor two for small Mach 
numbers, while for Mj > 5, the errors are smaller than 
30%. This could result in significant errors with the iden- 
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Fig. 3. Map of Lorentz factor distribution of a portion of the jet at a time before disruption, where a large amplitude 
wave is apparent (top panel, t = 320 Rj/c) and at the last frame (bottom panel, t — 1097 Rj/c). Coordinates are in 
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Fig. 4. Transversal structure of pressure perturbations (lighter shading indicates higher pressure values). Solid line 
indicates v z = 0.8 c contour. Three left panels: Cuts at 35 Rj, for t = 70, 140, 200 Rj/c where the elliptical or double 
helical mode rotation is apparent. Three right panels: Cuts at 105 Rj, for t = 210, 220, 240 Rj/c where the helical 
mode rotation is apparent. 



tification of the modes. Another difference we find comes 
from the fact that we do not observe the elliptical surface 
mode in the simulation, although it is fitted from observa- 
tions in LZ01 and we find in the solution to the stability 
problem that it has a short growth length (Fig. |5J). This 
could be due to the radial structure of the initial pertur- 
bation (oc sin(7rr)/ cosh(r)) giving zero initial amplitudes 
at the jet surfaces, and therefore suppressing the surface 
modes. Otherwise we would expect the elliptical surface 
mode to dominate at z < 60 Rj and maybe also farther 
downstream, as seen in the growth rates shown. Somehow, 
however, the helical surface mode seems to develop at dis- 
tances z > 70 Rj , and we think that this may be due to 
slight changes in the radial structure of the perturbations 
with downstream evolution as they grow in amplitude and 
modes interact among them. 

In the frame of the comparison between results from 
the simulation and linear analysis and from the fits in 
LZ01, we now focus our discussion on the three main struc- 
tures observed in the simulation. We define Af m = 4i?j, 
A™ = 25 Rj and Xf m — 50 Rj as the characteristic wave- 
lengths in the simulation. Propagation speeds of the per- 
turbations can be measured in the section of the jet dom- 
inated by the linear growth of instability. Although this is 
difficult due to the sparsity of the data frames, we derive 
the wave speed of the disruptive mode (A| im ) following 
the motion of the large amplitude wave (see Fig. |3J) from 
frame to frame, which gives u„,,2 ~ 0.38 c, or any fraction 



of this number. From Fig.[3wc can tell that the mode must 
be moving at v w ~ 0.09c. Another wave speed of the sys- 
tem is obtained from the rotation of the elliptical patterns 
similar to those shown in Fig. 01 which yields v w ~ 0.2 c. 
Unfortunately, we have not been able to find the wave- 
length of this elliptical mode in the pressure perturbation 
plots. We note that this speed is close to the value of v w 
given in LZ01 and it is different from the one predicted by 
the solution of the stability problem (Fig. 0) for the heli- 
cal first body modes that we claimed to generate this long 
scale structure as a beating pattern. This can be due to the 
presence of a thick shear layer that certainly changes the 
picture of the solution to the stability problem. It could 
also be due to the structure that we have used to mea- 
sure the speed having been artificially generated by the 
interaction of helical modes in 180° phase, and therefore 
giving a different velocity to those of single modes. Also, 
we use, as a limit for small wavelength perturbations, a 
wave speed equal to that of the flow (v Wt i ~ 0.88 c). 

The parameters of the simulated and resulting ob- 
served structures are given in Table [S] To reconcile the 
simulations with the observational results we calculate 

X obs from X sim (ytheor m Eq [gj using thc three d j ffer . 

ent values of v w mentioned in the previous paragraph. It 
should be noted that the two longest modes identified in 
the simulation have an observational counterpart when we 
take v w = 0.09c, which is a value close to that given by 
the stability problem for most of the identified modes in 
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Fig. 6. Theoretical computation of the pressure pertur- 
bation produced by the combination of modes identified 
in the simulation, compared to Fig. [SJ Lines indicate cuts 
at different distances from the axis. Dotted line stands 
for pressure perturbation at Rj/8, dashed line at 3Rj/8, 
dash-dot line at 5Rj/8, and long dashed line for 7 Rj/8. 
The wavelengths and wave speeds of the modes applied 
in this plot are indicated with arrows in Fig. [21 These are 
the helical first body mode at A = 4.53 Rj, A = 4.04 Rj 
for z < 60 Rj, the helical surface mode at A = 24.3 Rj 
and first body mode at A = 2.47 Rj and the elliptical 
first body mode at A = 3.22 Rj for z > 70 Rj. We have 
used a linearly growing amplitude between z 



70 Rj and 



120 Rj for the helical surface mode, between z 



70 Ri 



and z = 105 Rj for the helical first body mode, and con- 
stant amplitude from those distances. We have also used 
constant amplitude at z < 70 Rj for the elliptical body 
mode. 



\| OTi mode could be iden- 



the simulation (see Fig.|5J). The A 
tified with the fitted helical second body mode with the 
wavelength of 2 mas in LZ01. The Xf m mode could cor- 
respond to the fitted first body modes at the wavelength 
of 4 mas in LZ01. However, \™ m is identified as a helical 
surface mode, thus not coincident with the fitted second 
body mode to the 2 mas structure and we have used the 
long envelope of the beating structure in order to derive a 
4 mas structure, which is interpreted as two helical body 
modes with a similar wavelength in LZ01. Using the given 
wave speed of v w = 0.09c, the term sin^ has larger in- 
fluence on the result than the term 1 — v w cos 9j in Eq. 
so only larger angles to the line of sight would give larger 
X obs from a given X svm : the ARj wavelength would result 
in a 1 mas observed structure at 45° and the 25 Rj wave- 
length would result in a 6 mas observed structure at the 
same angle to the line of sight. Nevertheless, this angles 
are ruled out by jet to counter-jet (still not observed) flux 
ratios and by recent observations by Jorstad et al. ( 2005 ) . 
At that wave speed, shorter angles to the line of sight 
would result in even shorter observed wavelengths. This 



translates into the need of larger wavelengths in the simu- 
lation in order to fit them to observations, but Fig. [2]tclLs 
us that we are in the longest unstable wavelength limit for 
body modes, so this seems unrealistic. 

Why we do not see the 12 mas elliptical surface mode 
is thought to be due to the radial structure of the ini- 
tial perturbation (oc sin(7rr)/ cosh(r)) which, as stated 
above, gives zero initial amplitudes at the jet surfaces, 
therefore suppressing these modes. Moreover, the 12 mas 
mode, with the wave speed given in LZ01, would require a 
110 Rj wavelength in the simulation, which is difficult to 
observe even in a grid as large as was used here, in par- 
ticular when shorter harmonics grow fast and disrupt the 
flow. 



2.3.3. Nonlinear regime 

Nonlinear effects become important at time t — 350 Rj/c 
with the disruption of the head of the jet due to the longest 
helical mode A|"™ ~ 25 Rj . After that point perturbations 
produced at the disruption point propagate backwards 
slowly as a backfiow. The disruption point itself moves 
downstream due to constant injection of momentum at 
the inlet and the change of conditions around the jet. The 
disruption point advances from 160 Rj to 180 Rj by the 
end of the simulation (see Fig.[3J). 

Morphology of the jet at the end of the run (lower 
panel of Fig. |3J) is thus different from the observed source 
mainly due to the disruption of the jet. These difference 
may result from the development of a disruptive mode in 
the simulation which is not present in the real jet due 
to, for example, magnetic fields or an opening angle in 
the jet, not taken into account here. Uncertainties in the 
calculation of the physical parameters from the character- 
istic wavelengths following Eqs. Q and (0), as discussed 
in previous paragraphs, can be a source of error in the de- 
termination of the parameters of the jet, which, in turn, 
influence the long term stability properties. 

Disruption of the jet in the simulation contradicts ap- 
parently the fact that the jet is observed on much larger 
scales. It should be noted however that the disruption 
point still propagates outward at v ~ 0.03 c, at the end of 
the simulations. This implies that the simulation has not 
run long enough to reach a quasi-steady stage. It is also 
possible that the disruption observed is a transitory phase 
and that the simulation should have run longer in order to 
allow the jet to move downstream. If a stabilizing factor 
is needed in order to explain the jet in 3C273, we suggest 
several possibilities: 1) a thicker shear layer (Birkinshaw 
HTSM! Hardee & HuguesEEH Perucho et aLEEi, 2) in- 
clusion of the superluminal components in the simulations, 
as faster jets are much more stable against K-H instability 
(see Perucho et al. I2004bl Perucho et al. EHEJl, 3) a de- 
creasing density atmosphere (Hardee 119821 119871 Hardee 
et al. 120051 where it is shown, in the case of the jet in 
3C 120, that the expansion of the jet provides a stabiliz- 
ing influence), which must be the case as can be derived 
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Table 5. First column give identified wavelengths in LZ01 from shorter to longer, in the second column we have written 
characteristic wavelengths in the simulation also in increasing order, and the last three columns give the wavelengths 
as observed depending on the wave speed. 



X obs [mas] 


\ slm [Rj] 


c [mas] 


^v/%0.23 c [mas] 


c [mas] 


2 


4 


0.36 


0.44 


2.27 


4 


25 


2.28 


2.7 


14.3 


12 


50 


4.56 


5.5 


28.5 



from the outward dimming due to adiabatic expansion of 
the observed jet in the parsec scale , 4) a stabilizing con- 
figuration of magnetic field (Rosen et al. ^999 Frank et 
al. Jones et al. 115571 Ryu et al. 125551 Asada et al. 

2002 ) . The cumulative effect of these factors would effec- 
tively make the jet more stable already on the time scales 
probed by the simulation. 

2.4. Simulation 3C273-B 
2.4.1. Initial setup 

In this simulation, the stationary flow is perturbed by pre- 
cession and periodical ejections of faster components. The 
initial conditions are similar to those in the first simula- 
tion. The width of the shear layer is ~ 0.6 Rj (m = 8). 
The precession frequency is derived from the observed 
~ 15 yr (6 Rj / c in the code units) periodicity of posi- 
tion angle variations (Abraham et al. 11996(1 . The fre- 
quency of ejections of components is set by the reported 
~ lyr (OA Rj/c) periodicity in the optical light curve 
(Babadzhanyants and Belokon 1993). The duration of each 
ejection is estimated to be ~ 2 months (0.066 Rj/c), set by 
the approximate inspiralling time from an orbit at w 6i?c 
around a 5.5 x 10 8 M Q black hole. The amplitude of the 
precession is set by the true opening angle (0.4°) of the 
jet obtained by deprojecting the apparent opening of 1.5°. 
This precession is included in transversal velocities at the 
injection point as follows: 

v± = A v z (cos(ujt), sm(ut)) (7) 

where v± and v z are the transversal and axial components 
of the velocity, t is time, Aq = 6.83 10 -3 is the initial 
amplitude, and w ~ 1.01 c/Rj is the angular frequency 
calculated from the precession period of ~ 15yrs. 

The fluid in the injected components is considered to 
have the same density as the underlying flow and to be 
in pressure equilibrium with it. Velocity of the fluid in 
components is taken as constant, with the mean value 
of Lorentz factor j c ~ 5 as reported in Abraham et al. 
(1995J). The components are generated as shells of fluid 
with a diameter of 0.5 Rj ejected along the axis. 

The numerical grid for this simulation covers 30 Rj 
(axial) times 6Rj times 6 Rj (transversal), i.e., 24 pc x 
4.8 pcx 4.8 pc. The axial dimension of the grid is related to 
the ejection of the components. We take into account that 
the wavelength induced by the precession of components, 
if they move ballistically, is A c ~ Pv c , where P is the 



precession period and v c is the injection velocity of the 
fluid in the components. This gives A c ~ 6Rj, and we 
have chosen the grid of 5 A c to allow the wave to become 
apparent. 

The resolution of the grid is 16 cells/ Rj in the transver- 
sal direction and 32 cells/ Rj in the direction of the flow. 
An extended grid is introduced in both transversal and 
axial directions. In the radial directions, it has 36 cells 
reaching out to 15 Rj on each side of the jet (increasing 
the cell size by 7.7% from one cell to the next). In the ax- 
ial direction the extended grid has 192 cells, reaching up 
to 45 Rj . This simulation has lasted for a time of 70 Rj / c 
(i.e., more than two light crossing times of the grid). 

2.4.2. Results 

Fig. [7J shows the solutions of the linear stability prob- 
lem for the jet parameters given in Tableland indicates 
the characteristic wavelengths arising un this simulation. 
Fig. [S] shows longitudinal cuts of pressure perturbation at 
different jet radii at time 40 Rj/c. Inspection of the lon- 
gitudinal cuts indicates that, close to the injection point 
and to the jet axis, a symmetric, short wavelength pertur- 
bation generated by the fast components dominates the 
structure of the flow. Its wavelength is Xf m — 0.4 Rj, 
and it is clearly related to the ejection period of com- 
ponents (0.4 Rj/c). After 3 — 4i?j, the presence of the 
fast components is also evident in the jet boundary (at 
r = 7/&Rj). These high-frequency and short-wavelength 
structures are damped at z > 10 Rj, as expected to oc- 
cur for perturbations with wavelengths smaller than the 
shear layer width. Close to the jet boundary, the most pro- 
nounced structure is the typical antisymmetric pattern of 
helical motion induced by the precession (A|™ ~ 3.7 Rj). 
This structure is driven by the precession of the injected 
components and it couples to the helical surface mode at 
a frequency of 0.8c/ Rj. A pinch mode structure with a 
wavelength Xf" 1 = 3.7 Rj is also observed at r = Rj/8, 
coinciding with the maximum growth of the first body 
mode (see Fig.|7Jand Table EJ). 

Farther downstream and close to the jet center, there 
is a Xf m = 0.9 Rj wavelength on top of the longer pinch 
mode (Ajf m ). This short wavelength coincides with the 
maximum growth rate of the elliptical third body mode 
at the helical driving frequency (1.01 c/Rj) (see Fig. [7J 
and Table 0}. We remind the reader that this is frequency 
correspond to a 360° turn, as explained in the caption 
of Fig. The radial structure of this mode in the simu- 
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Fig. 7. Solution of the stability problem for the parameters given in Table [21 with the characteristic frequencies 
for simulation 3C273-B. Frequency ujh = 1.01c/ Rj corresponds to that excited in the simulation 3C273-B. Arrows 
connecting the wavelength and wave speed plots indicate identified modes in the simulation 3C273-B (see text). 



lation is also coincident with the theoretical structure of 
the third body elliptical mode, with the maximum am- 
plitude at r — Rj/4: and decreasing amplitude at larger 
radii. The helical driving frequency falls very close to the 
maximum growth of the second body helical mode with 
wavelength Af"™ = 1.5 Rj, observed in Fig. |5] However, 
the radial structure found in the simulation is somehow 
a mixture of the second body mode with that of the first 
body mode. We expect the second body mode to have 
an amplitude maximum at r = 1/SRj, but we find that 
this maximum occurs at r = 3/8Rj. The first body mode 
at this wavelength develops an amplitude maximum at 
r = 3/8Rj, which is coincident with the one found in the 
simulation, although the amplitude at r = 1/SRj found 
in the simulation is too large for this mode. We interpret 
this as the second body mode being triggered at the driv- 
ing frequency, which in turn excites the first body mode 
at the same wavelength. Both modes seem to be triggered 
out of phase and interfere destructively in the inner jet, 
but the first body mode dominates in the mid jet region 
and both decline in amplitude towards the jet surface. 
Comparison of smaller and larger radius plots of pressure 
perturbation in Fig.[H]shows large positive offsets observed 
at larger radii. This indicates possible drifting of compo- 
nents (shocks) to outer radii of the jet as they follow the 
helical path given by the surface mode (similar to what 
has been reported by LZ01). This feature would produce 
enhanced emission regions at the positions of the helix 
in which the flow moves in a direction closer to the line 



of sight. Using all of the modes identified in the simula- 
tion we have computed a theoretical representation of the 
pressure perturbation in the jet. Fig.UHshows the resulting 
axial cuts of the pressure perturbation and indicates that 
the theoretical calculation generates most of the structures 
found in the simulation, except those which are intrinsi- 
cally nonlinear (for instance, the injected components at 
the jet inlet and those structures at outer jet radii farther 
downstream) . 

The complexity of the structure is further illustrated 
by the surface plot of the flow Lorentz factor shown in 
Fig. ^3 Although in this figure we only see patterns in 
the fluid, comparison of the wavelengths seen here and 
in Fig. allows us to identify the wavelengths derived 
from the fluid patterns in Fig. 1101 and from the frequency 
of injection of components. The middle panel of Fig. 
(Lorentz factor 7 = 2.5) indicates that the periodicity in- 
duced by individual jet components dominates the struc- 
ture up to distances z ~ 10 — 15 Rj, but farther down- 
stream the components expand longitudinally and start 
to interact with each other, generating a semi-continuous 
structure that is dominated by the helical motion induced 
by the precession. The top panle of Fig. ^3 (Lorentz fac- 
tor 7 = 3.0) indicates that the distinct regions of the flow 
moving at higher speed disappear downstream. This can 
be explained by the deceleration of the fluid in the compo- 
nents. The deceleration can be caused either by the inter- 
action with the background flow, or by a radial and lon- 
gitudinal expansion. These results are in agreement with 
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Fig. 8. Longitudinal cuts of pressure perturbation at dif- 
ferent radii (Rj/8 top left, iRj/8 top right, 5Rj/8 bottom 
left, 7Rj/8 bottom right) and in symmetric positions with 
respect to the jet axis at t = 40Rj/c. Short wavelength 
(Af m = 0.4 Rj and Xf m = 0.9 Rj) structures are observed 
close to the jet axis. The amplitude of the first structures 
decrease downstream as a consequence of the drift of com- 
ponents to outer radial positions due to the precessional 
motion and their interaction with the underlying flow. 
Also at positions close to the axis, we find a structure 
with wavelength Xf m = 1.5 Rj. Short wavelength sym- 
metric (Af lm = 0.4 Rj) and longer wavelength antisym- 
metric (A;f m = 4Rj) structures are also observed at outer 
radial positions. The longer structure dominates close to 
the injection, before the components expand and/or drift 
from the axis to the jet boundaries, when short scale struc- 
tures can be observed to be modulated by the longer one 
(from z <~ 4:Rj). After this, at z — 9Rj, the antisymmet- 
ric longer structure grows in amplitude, with some spikes 
that can be associated with the drift of the components 
due to their precessional motion. 



results by Lobanov and Zensus l|1999jl . and Lobanov and 
Roland (2001), suggesting that shocks dominate the jet 
structure close to the nucleus, and that fluid instabili- 
ties become important farther downstream. On the other 
hand, it is not clear at the moment whether stronger fluid 
components (i.e., faster and denser) could survive longer 
in the jet. We finally note that the precession wavelength 
obtained from the simulations is smaller than the one cal- 
culated theoretically from the advance speed of compo- 
nents (Ajf m = ARj versus A c ~ QRj), which maybe taken 
as an indication of non-ballistic motion of the components. 

In TableElwe list possible observed wavelengths corre- 
sponding to the two main wavelengths identified from the 
simulation. It is clear that the observed wavelength of pre- 
cession (18 mas) cannot be recovered even with extremely 
fast components (v w — c gives \ obs = lOmas) with 
the adopted viewing angle of 15°. Moreover, if we con- 
sider the mean apparent proper motion of 0.93 /j _1 mas/yr 
(Abraham et al. I199fl[l . the apparent speed is f3 app ~ 10. 
This speed cannot be reconciled with a viewing angle of 
15°, as the resulting intrinsic speed is larger than 1, and it 
would require a Lorentz factor of 7 = 10 if the viewing an- 




Fig. 9. Theoretical computation of the pressure pertur- 
bations produced by the combination of modes identified 
in the simulation. Different panels indicate cuts at differ- 
ent distances from the axis in the X-Z plane. The wave- 
lengths and wave speeds of the modes contributing to this 
perturbations are marked by arrows in Fig. \7\ These are 
the first pinch body mode at A — 3.7 Rj, the helical sur- 
face mode at A = 3.7 Rj, helical first and second body 
modes at A — 1.5 Rj, and the elliptical third body mode at 
A = 0.9 Rj. We have used a linearly growing amplitude for 
the modes. This perturbations can be compared with the 
structures identified in Fig. |H1 Main differences between 
this plot and Fig. [S] are due to nonlinearities introduced 
by injected components, mainly close to the injection on 
the jet axis and at outer radii farther downstream. 



gle is reduced to 10°. Recent work by Jorstad et al. (2005) 
indicates that the jet in 3C273 can have a viewing angle 
as small as 6° and component Lorentz factors of 7 ~ 10.6. 
These numbers would transform A|"™ = 4 Rj into an ob- 
served wavelength 14.1 mas. This wavelength is in fair 
agreement with the 18 mas mode assigned to precession 
in LZ01. In summary, either the precession period should 
be longer, or the viewing angle should be smaller than 
15° and component Lorentz factors higher than 7 = 3, in 
order to reconcile the 18 mas structure with precession. 

Another question addressed by this simulation is 
whether the periodic injection of fast components could 
generate smaller structures observed by LZ01 (the 2 mas 
and 4 mas modes) , where these wavelengths are identified 
with the elliptical modes of Kelvin-Hclmholtz instability. 
In our simulation, the fast components generate mainly 
pinching modes, although this is simply due to their sym- 
metric nature. Again, we find that the structures gener- 
ated in the simulation are small compared to those ob- 
served. However, in the simulation, the fluid in the com- 
ponents moves with a smaller speed than the injection 
one (7 = 5, see Fig. I10|l and therefore, smaller than that 
found from the observations (7 ~ 5 — 10). We relate this 
fact to the possible slowing of the components themselves 
as they propagate downstream. As pointed out before, a 
possible cause for this could be that the simulated fluid 
injected in components has the same density as the back- 
ground flow, whereas if it was denser (as could be expected 
as generated from strong accretion activity) or propagat- 
ing in a decreasing density atmosphere, this fluid would 
have larger inertia and would generate faster components. 
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Table 6. First column gives the observed wavelengths in the simulation, and the last three columns give this wave- 
lengths as observed depending on the wave speed. We have used, for the wave speed, some of those obtained from 
Fig.0for the identified modes (second to fourth columns), the underlying flow speed (fifth column) and the maximum 
speed from the components (sixth column). 
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Fig. 10. Surface plots for three different values of Lorentz factor: Upper panel, 7 = 3.0; middle panel, 7 = 2.5; lower 
panel, 7 = 2.0. The precession motion, coupled to a helical mode, is more apparent on the slower, underlying flow, 
moving at 7 = 2.0 — 2.5, as shown in LZ01. The surface plot for 7 = 3.0 illustrates the appearance of the jet at higher 
frequencies and smaller viewing angles, where the emission is dominated by a more relativistic plasma. 



We find that, in order to produce a 4 mas structure, we 
need the Lorentz factor of components to be 7 ~ 30, 
whereas 7 ~ 15 is required to explain the 2 mas wave- 
length, if we keep the 1 yr period. Adopting the longest 
measured period of the ejections of 1.7 yrs (Abraham and 
Romero 1999J , these values would be reduced to 7 ~ 17 
and 7 ~ 8.5, respectively. The same authors gave a pe- 
riodicity in the injection Lorentz factor of about 4 yrs; if 
we consider this period as the generator of short modes, 
7 ~ 7.5 and 7 ~ 4 could explain those structures. The lat- 
ter values agree well with the Lorentz factors inferred from 
the observed kinematics of the jet. This means that the 2 
and 4 mas wavelengths should be associated only with the 
strongest and fastest ejections occurring roughly once ev- 



ery 4 yrs. With the jet parameters given by Jorstad et al. 
( 2005 ), the inferred wavelength for the shortest symmet- 
ric structure (OARj) is 1.4 mas, which is within a factor 
of 3 from the 4 mas wavelength identified in LZ01 with 
the elliptical first body mode. The symmetric or pinching 
nature of the perturbations induced by such components 
could not explain the double helix structure in the jet, but 
we observe in the simulation that elliptical modes can be 
triggered by the presence of pinching and helical pertur- 
bations. 

3. Conclusions 

We have performed two numerical RHD simulations with 
different initial setups in order to study the physical pro- 
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cesses generating the observed structures in the parsec- 
scale radio jet in the quasar 3C 273. In the simulation 
3C273-A, we have included a general set of helical and 
elliptical perturbations in a long jet with the basic phys- 
ical parameters adopted from LZ01. In the simulation 
3C273-B, we have used a shorter jet with the same physi- 
cal parameters and have included precession and injection 
of fast components. The simulation 3C273-A was aimed 
to generate structures with wavelengths similar to those 
measured by LZ01 from the growth of Kelvin-Helmholtz 
perturbations. The simulation 3C273-B was designed to 
check if by combining the ejection of superluminal com- 
ponents and jet precession, with the periodicities reported 
in Babadzhanyants and Belokon (1993) and Abraham et 
al. l|199rj[) . the same structures could be generated. 

We find that the structures generated in simulation 
3C273-A are of the same order in size as those observed, 
if the relativistic propagation effects of the waves are taken 
into account. We observe in the solution of the stability 
problem that the instability modes found in the simulation 
propagate at mildly relativistic speeds. These wave speeds 
differ from those derived from the approximations used in 
LZ01. This can be due to the uncertainties introduced by 
the approximations to the characteristic wavelengths in 
the interpretation of the observations in LZ01. However, 
we show that wavelengths similar to the observed ones 
are found for the wave speed given by the solution of the 
linear problem, although the modes fitted in LZ01 and 
those used here for the same wavelengths are not coinci- 
dent. The solutions of the stability problem applied to the 
adopted wave speeds (0.23 c) and line of sight (15°) show 
that any body modes present in the jet should be much 
shorter than those fitted in LZ01. It should be noted that 
these differences do not rule out the presence of Kelvin- 
Hclmholtz instability in parsec-scale jets. Despite difficul- 
ties in the mode identifications, the structures generated 
in the simulation are similar to those observed by LZ01. 

Regarding the long-term stability of the flow, we note 
that the jet in the simulation 3C273-A is disrupted at 
~ 170 pc from the inlet, contrary to the observations trac- 
ing the jet in 3C273 up to 60 kpc away from the source. 
The reasons for this difference may be found in the con- 
junction of several factors. 1) The numerical simulation 
does not run long enough to reach a fully steady-state 
regime. The disruption point moves downstream along 
the simulation, which could imply that the disruption is a 
transitory phase. 2) Magnetic fields have not been taken 
into account neither in the linear analysis, nor in the nu- 
merical simulation - and it should be noted that the mag- 
netic fields may be dynamically important at parsec scales 
(Rosen et al. Frank et al. HTMl Jones et al. ITM71 

Ryu et al. I2UUU1 Asada et al. I2UU2|) . 3) We only simu- 
late the underlying flow, without considering the faster 
and possibly denser fluid in the superluminal components. 

4) Inaccuracies in the linear analysis approximations can 
lead to large uncertainties in physical parameters derived. 

5) Differential rotation of the jet, shear layer thickness 
CBirkinshaw UMTl Hardee & Hueues EHoU Perucho et al. 



2006$, and a decreasing density external medium could 
also play an important role (implying jet expansion; see 
Hardee HnH2IIM3 and Hardee et al.EDDSJ. 6) Arbitrary 
initial amplitudes of perturbations were chosen for the 
simulation, so we could have included too large pertur- 
bations. A combination of these factors could well change 
the picture of the evolution of the jet in terms of its sta- 
bility properties. The effects of the rotation and magnetic 
fields on the stability of jets remain unclear, since no sys- 
tematic numerical study has been performed up to now. 

In the simulation 3C273-B, we studied the effect of 
precession on the jet evolution and investigated the possi- 
bility that the short wavelength structures found in LZ01 
were not due to K-H instabilities but due to the period- 
icities induced in the flow by the ejection of components. 
We demonstrate that such non-linear features as superlu- 
minal components generate linear structures in the form 
of Kelvin-Helmholtz instabilities which can be analyzed 
in the framework of linear perturbation analysis. One of 
the main conclusions that can be derived from this sim- 
ulation is that helical twists can be excited by periodic 
injections if there is some induced hclicity in the system. 
This hclicity is induced in our simulation by the helical 
perturbation frequency, but in real jets this hclicity could 
be induced by helical jet magnetic fields and/or by jet 
rotation. We have also shown that, in order to explain 
the observed 18 mas wavelength in terms of precession, ei- 
ther longer driving periodicities than the 15 yr suggested 
by Abraham & Romero (1999) would be needed, or this 
wavelength must be induced by very fast components ob- 
served in a jet moving at a viewing angle 8 < 15°. The fast 
components could only generate the shorter wavelengths 
given in LZ01 (A = 2 mas and A = 4 mas) if a proper 
combination of the velocities and injection periodicities is 
used. Altogether, we find that inclusion of faster compo- 
nents and precession with the 15 yr periodicity does not 
explain well the observed wavelengths and periodicities. 
This gives more weight to the general conclusion about 
K-H instability acting prominently in the flow. 

In the future, numerical simulations of this kind may 
be used to constrain the basic parameters of the flow such 
as the viewing angle and the component speed. Inclusion 
of magnetic fields, differential rotation and the effects of an 
atmosphere with a decreasing density could help reconcil- 
ing better the simulations with the observed structures. 
In this way, for example, an increase of the jet radius 
due to decreasing external density could cause a down- 
stream increase of wavelengths of K-H instability modes 
(Hardee et al. l2005[) . This remains to be seen with future, 
full-fledged RMHD simulations of the relativistic jet in 
3C 273. The scope of the present work could be expanded 
to other sources, and applied to prominent jets for which 
the transversal structure may be resolved, such as 3C 120, 
extending the work done by Hardee et al. I|2005|) by per- 
forming numerical simulations. 
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